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ABSTRACT 

T Cha is a nearby (d ~ 100 pc) transition disk known to have an optically thin gap separating optically 
thick inner and outer disk components. Huelamo et al. (2011) recently reported the presence of a low- 
mass object candidate within the gap of the T Cha disk, giving credence to the suspected planetary 
origin of this gap. Here we present the Herschel photometry (70, 160, 250, 350, and 500 /im) of 
T Cha from the "Dust, Ice, and Gas in Time" (DIGIT) Key Program, which bridges the wavelength 
range between existing Spitzer and millimeter data and provide important constraints on the outer 
disk properties of this extraordinary system. We model the entire optical to millimeter wavelength 
spectral energy distribution (SED) of T Cha (19 data points between 0.36 and 3300 fim without any 
major gaps in wavelength coverage). T Cha shows a steep spectral slope in the far-IR, which we find 
clearly favors models with outer disks containing little or no dust beyond ^^40 AU. The full SED can 
be modeled equally well with either an outer disk that is very compact (only a few AU wide) or a 
much larger one that has a very steep surface density profile. That is, T Cha's outer disk seems to 
be either very small or very tenuous. Both scenarios suggest a highly unusual outer disk and have 
important but different implications for the nature of T Cha. Spatially resolved images are needed to 
distinguish between the two scenarios. 

Subject headings: circumstellar matter — protoplanetary disks — submillimeter: planetary systems 
— planetary systems: planet-disk interactions — stars: individual (T Cha) 



1. INTRODUCTION 

Transition disks were first identified by the Infrared 
Astronomical Satellite as objects with little or no excess 
emission at A < 10 /im and a significant excess at A > 10 
//m (Strom et al. 1989). This broad definition encom- 
passes a wide range of SED morphologies (Najita et al. 
2007). T Cha belongs to a subclass of transition disks 
known as cold disks, which are characterized by a steep 
rise in the SED between ~10 and ~ 30 jam (Brown et al. 
2007; Merm et al. 2010). It is now well established, from 
both SED modeling and direct (sub) millimeter imaging, 
that the peculiar SED shapes of cold disks are due to the 
presence of large inner holes or gaps in their dusty disks. 
Inner holes can be produced by a number of processes, 
including grain growth, photoevaporation, and dynami- 
cal interactions with (sub) stellar or planetary mass com- 
panions (see Williams & Cieza, 2011 for a recent review). 
However, objects such as T Cha, LkCa 15, and UX Tau 
are believed to have an optically thick inner disk within 
a much larger inner hole (Olofsson et al. 2011, Espail- 
lat et al. 2010). Optically thin gaps separating optically 
thick inner and outer disk components are diflicult to 
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explain by anything other than the dynamical interac- 
tion of substellar or planetary mass objects embedded 
within the disk. The discovery of what seems to be ei- 
ther a young planet or a brown dwarf surrounded by a 
significant amount of dust within the gap of the T Cha 
disk (Huelamo et al. 2011) supports the suspected plan- 
etary origin of the gap. Detailed studies of the properties 
of the T Cha disk thus become highly desirable to fur- 
ther our understanding of protoplanetary disk evolution 
and planet formation. Using Spitzer photometry and 
spectroscopy. Brown et al. (2007) first modeled T Cha 
with an inner disk extending from 0.08 to 0.2 AU and an 
outer disk extending from 15 to 300 AU. More recently, 
Olofsson et al. (2011) simultaneously modeled the SED 
of T Cha and near-IR interferometric observations ob- 
tained with the AMBER instrument at the Very Large 
Telescope Interferometer (VLTI; Scholler 2007). They 
spatially resolved the narrow inner dusty disk and con- 
strained its inner edge to be ~ 0.1 AU from the star. 

Here we present the Herschel photometry (at 70, 160, 
250, 350, and 500 /im) of T Cha and model its entire op- 
tical to millimeter wavelength SED. We focus our anal- 
ysis on the properties of the outer disk (i.e., beyond the 
gap), which is the region that is both the best probed 
by Herschel data and the least constrained by previous 
studies. 

2. OBSERVATIONS 

2.1. Herschel observations and data reduction 

We obtained far-IR and submillimeter wavelength pho- 
tometry for T Cha as part of the DIGIT Open Time 
Herschel Key Program using both the Photodetector Ar- 
ray Camera and Spectrometer (PACS, 70 and 160 /xm, 
Poglitsch et al. 2010) and the Spectral and Photometric 
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Imaging Receiver (SPIRE, 250, 350 and 500 /im. Grif- 
fin et al. 2010). The PACS observations (OBS ID = 
1342209063) were obtained on November 5^^, 2010, us- 
ing the Mini Scan Map mode, a repetition factor of 3, 
and a scan angle of 70 degrees. Similarly, the SPIRE ob- 
servations (OBS ID = 1342203636) were taken on August 
24^^, 2010, using the Small Scan Map mode, and a repe- 
tition factor of 3. Both the PACS and SPIRE data were 
processed using HIPE (Herschel Interactive Processing 
Environment; Ott et al. 2010) version 7.1. The PACS 
beam sizes at 70, and 160 /im are 5.5^' and 11^', and the 
data were resampled to and 2'^ pixel images, respec- 
tively. SPIRE beam sizes at 250, 350 and 500 /im are 
18.1, 25.2 and 36. 6'^ respectively, and the images were 
resampled to 6'^ 10'^ and 14^' pixel images, respectively. 

T Cha was detected with good signal to noise ( > 15- 
50) at all PACS and SPIRE wavelengths. PACS and 
SPIRE fluxes were derived with a psf-fltting routine es- 
sentially identical to the c2dphot software used for the 
c2d Spitzer Legacy Program (Evans et al. 2007). We 
used the empirical psf's posted by the NASA Herschel 
Science Center'^ and found that they gave a good match 
to the actual images of T Cha. In addition, we checked 
that aperture photometry gave similar results and as- 
signed the flnal uncertainties based on the 10% and 15% 
error estimated for PACS and SPIRE absolute photome- 
try, respectively. We also added an additional 5% uncer- 
tainty to reflect incomplete knowledge of the psf and the 
agreement between the psf-fltting and the aperture pho- 
tometry. The resulting Herschel photometry for T Cha 
is listed in Table 1. 

2.2. Photometry data from the literature 

In order to construct the full SED of T Cha, we col- 
lected optical, near-IR, mid-IR, and mm wavelength pho- 
tometry from the literature, which are also listed in Table 
1. The full optical to mm wavelength SED for T Cha is 
shown in Figure [T] There are now a total of 19 photom- 
etry points without major gaps in wavelength coverage 
between 0.36 and 3300 /im. In addition to the above pho- 
tometry, the T Cha SED is well sampled by a Spitzer-IRS 
spectrum, covering the 5.2 to 38 /im region (Brown et al. 
2007). The sharp dip in the SED around 15 /im indicates 
the presence of a wide gap in the disk. 

3. DISK MODEL 

The starting point of our study is the model presented 
by Olofsson et al. (2011, Oil hereafter). Following Oil, 
we use the Monte Carlo radiative transfer code MCFOST 
(Pinte et al. 2006; 2009). We parameterize the struc- 
ture of the T Cha disk with the following parameters: 
the inner and outer radii (R^^ and Rout^ respectively), 
the index a for the surface density proflle (S (r) = S 
(r/ro)" ) and a disk scale height H{r) assuming a verti- 
cal gaussian distribution (exp[— z^/2H(r)^]), appropriate 
for a disk in hydrostatic equilibrium. The disk's flaring 
is described by a power-law determining the scale height 
as a function of the radius {H{r) = Ho(r/ro)^). The 
dust content is described by a differential power-law for 
the grain size distribution (dn(a) ex a ^da), between the 
minimum {^rnin) and maximum (a 

max 

) grain sizes. The 

total mass of dust in the disk with grain sizes between 

^ https:/ /nhscsci. ipac.caltech.edu/sc/index.php/ 



8imin and 8imax IS givcu by the Mdust parameter. Each 
of the above parameters can be deflned independently 
for the inner and outer disk components, which are sep- 
arated by a wide gap. 

The main goal of this Letter is to investigate the con- 
straints the Herschel data provide on the outer disk prop- 
erties of T Cha. As the Herschel data give no new in- 
formation on the inner regions, we adopt the inner disk 
parameters from Oil, which are listed in Table 2, as 
flxed values. The stellar parameters for T Cha are also 
taken from Oil and are listed in Table 2 as well. We 
also flx the minimum grain size for the outer disk and 
its inclination. We are thus left with 7 free model pa- 
rameters, all of which are related to the outer disk: R^^, 
Rout, Oi, ^max, ^^dust, and H50 (the scale height at 
50 AU). Motivated by the lack of 10 /im silicate emis- 
sion in the Spitzer-IR^ spectra of T Cha, Oil adopted 
for the inner disk a mixture of astro silicates (Draine 
& Lee 1984) and amorphous carbon (Zubko al. 1996) 
with a mass ratio of 4:1 between silicate and carbon. For 
simplicity, we adopt the same silicate and carbon mix 
for the dust composition of the outer disk, but remind 
the reader that disk parameters derived from SED fltting 
are strongly dependent on the choice of dust properties. 
Although discussing the dependence of each parameter 
on dust composition is beyond the scope of this Letter, 
we have verifled that our main conclusions remain un- 
changed if pure astro silicate grains are adopted for the 
outer disk. 

3.1. Model grid and best-fit models 

In order to constrain the parameter space consistent 
with the observed T Cha SED, we ran a grid of models 
varying the 7 free parameters in our model setup. The 
values sampled by the grid are listed in Table 3. The 
parameters j3 and H50 are sampled linearly, while R^^, 
and Rout, ol, ^max, and Mdust are sampled logarithmi- 
cally. The number of values per parameter range from 3 
to 10 and the total number of models amount to 126000. 
The parameters of the best-flt model (i.e., with the min- 
imum x^) are also listed in Table 3. The global best-flt 
model so derived for the outer disk of T Cha is a nar- 
row ring with an inner radius of ~18 AU and a width 
of only ~2 AU. This narrow ring remains optically thick 
even at (sub) millimeter wavelengths, which explains the 
very large mass that is allowed by the model. Its SED is 
shown as a solid black line in Fig. 1. 

Since SED modeling is known to be highly degenerate, 
the best-flt model is unlikely to be a unique solution. We 
therefore adopt the Bayesian method (Press et al. 1992; 
Pinte et al. 2008) to calculate the probability of different 
values for each parameter given the available data. We 
assume that we have no prior knowledge on the model 
parameter, which implies that the relative probability of 
a given model is proportional to exp[— x^/2]. All proba- 
bilities are normalized so the sum of all the probabilities 
of the models in the grid is equal to 1. The top two panels 
in Fig. 2 show the probability that each parameter takes 
a given value from the grid, calculated by summing over 
the probabilities of all the other parameters. Clearly, the 
^dust, ^max, and Rin parameters are much better con- 
strained than the others. The lower panels in Fig. 2 show 
the joint probabilities of the Rin, Rout, and a parame- 
ters. Models with R^^ ~ 18 AU and slightly larger Rout 
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values are clearly favored to best reproduce our observa- 
tions. There are a large number of narrow ring models 
that can reproduce the observed T Cha SED relatively 
well. The best-fit model is an instance of a large family 
of such models. This can be understood considering that 
the SED of a narrow ring is little sensitive to parameters 
such as a, /3, or even Mdust and ^raax- Wider rings are 
also allowed, but they require fine-tuning of other param- 
eters (see for example the red dashed line in in Fig. 1 and 
the ce = — 1 model in Table 3) and thus appear as less 
likely models in the joint probability plots. Large disks 
(Ront ^150 AU) with shallow surface density profiles 
{a > -1) result in a poor fit to the SED. This can be seen 
in Fig. 2 in the dark region in the upper right corner of 
the OL vs Ro^xt plot (lower row) and in the red histogram 
in the a probability distribution (middle row). We find 
that all such models result in far-IR Herschel colors that 
are much redder than observed (see bottom panel in Fig. 
2). In other words, all disk models with a shallow sur- 
face density profile {a > —1) and a large outer radius 
(e.g., > 150-300 AU) either under-predict the observed 
30 and 70 /im fiuxes and/or over-predict the 250 and 350 
/im fiuxes. This is not surprising considering that a nar- 
row ring around T Cha (like the best-fit model) would 
be "missing" most of the cool material emitting at 250 
and 350 fim expected for a large disk with a relatively 
shallow surface density profile. However, a narrow ring is 
not the only way to reproduce the far-IR Herschel colors. 
We also find that the outer disk can be arbitrarily large 
(e.g., Rout = 300 AU) as long as the surface density pro- 
file is sufficiently steep (e.g., a = —3) to quickly reduce 
the amount of dust at large radii. See for example the 
alternative a = —3 model in Table 3, corresponding to 
the light blue dot-dashed line in Fig. 1. This alternative 
model falls outside of our parameter grid, but illustrates 
well the strong degeneracy between the outer radius and 
the surface density profile of the disk. We note that while 
the a = —1 and a = —3 models have very different Rout 
values, their Mdust, Rm, a^ax, and H50 parameters 
are much more similar. 

Independently of the probability space they occupy in 
Fig. 2, the 3 models discussed above and listed in Table 
3 are fully consistent with all the available data and, in- 
dividually, can be considered to be equally good descrip- 
tions of the T Cha disk. One feature that our models 
cannot reproduce well is the pronounced V-shape of the 
Spitzer-IRS spectra around 15 /im. The models also un- 
derestimate the 2.2 and 3.6 /im excesses from the inner 
disk, which parameters are heavily driven by near-IR in- 
ter ferometry observations. Near-IR variability makes it 
difficult to simultaneously fit multi-epoch observations in 
this spectral region (see Oil for a detailed description of 
the inner disk modeling). Because the inner disk is op- 
tically thick, the exact geometry of the inner disk might 
affect some of the outer disk parameters like H50 and 
R^^ (Mulders et al. 2010), but is unlikely to play any 
role on the values of a and Rout, the key parameters in 
this study. 

4. DISCUSSION 

The results from the previous section show that T Cha 
SED can be modeled with multiple combinations of pa- 
rameter values. The most likely models can be divided 
into two broad categories: narrow rings and wider disks 



with steep surface density profiles. Both families of mod- 
els suggest a highly unusual outer disk and have impor- 
tant but different implications for the nature of T Cha. A 
narrow ring would strongly suggest that the outer disk 
has been truncated by the dynamical interaction of a 
companion, as is known to be the case in the HD 98800B 
and Hen 3-600 A multiple systems, which have outer radii 
of only 10 and 20 AU, respectively (Andrews et al. 2010). 
Very deep imaging of T Cha (Chauvin et al. 2010; Vi- 
cente et al. 2011) rule out any stellar companion that 
could be responsible for the truncation of the outer disk, 
unless its projected separation is currently less than 0.1^' 
(10 AU). Furthermore, recent aperture masking obser- 
vations by Huelamo et al. (2011) only report a very 
low-mass companion candidate (either a young planet 
or a brown dwarf) with a projected separation of 6.2zb0.7 
AU. Since this object is the most likely cause of the inner 
cavity (R^^ ~ 10-20 AU based on the SED modeling), its 
physical separation must be < R^n- Any object truncat- 
ing the outer disk thus must be substellar or planetary 
in nature. 

We speculate that the recently reported companion to 
T Cha could potentially also provide the explanation for 
the models with a < — 2. These unusual a values would 
imply that the disk material is highly concentrated to- 
ward the inner edge of the outer disk (i.e., that material 
is accumulating at the outer edge of the gap, at ~15 AU). 
As a reference, a = —1.5 in the canonical Minimum Mass 
Solar Nebula (Weidenschilling, 1977). T Cha's outer disk 
is relatively massive but shows little or no accretion onto 
the star (Alcala et al. 1993, Schisano et al. 2009). The 
transport of angular momentum and material across the 
outermost part of the disk should not be affected by the 
presence of the low-mass object, but when circumstel- 
lar material reaches the outer edge of the gap, it can be 
accreted by the low-mass object, channeled toward the 
inner disk, or accumulated at the gap's edge (Lubow & 
D'Angelo, 2006). The dynamical interaction between a 
gas-rich disk and a planetary mass object is a complex 
hydrodynamic problem that depends on the mass and lo- 
cation of the planet as well as on the properties of the disk 
(viscosity, disk scale height, etc), but an a value of < —2 
would suggest that it could have an important effect on 
the surface density of the outer disk. Recent 3-D hy- 
drodynamic models of planet-disk interactions show the 
accumulation of mm-sized grains in the pressure maxima 
of the gas caused by the spiral density waves triggered by 
the planet (Fouchet et al. 2010). This implies that the 
density enhancement at the edge of the gap of the large 
grains traced by Herschel could be even larger than it is 
for the gas. However, given the degeneracy between a 
and Ront, the dynamical effect of the companion on the 
disk cannot be tested without spatially resolved observa- 
tions. Our models show that, given the SED data alone 
and without knowing Rout, the value of a for the dust 
disk could be anywhere between —0.75 and —3. 

5. CONCLUSIONS 

New Herschel observations of T Cha reveal a steep 
spectral slope between 70 and 350 /im (i.e., very blue 
far-IR Herschel colors compared to a "normal disk"). 
The steeply decreasing far-IR SED is not due to a lack 
of large grains. In fact, the (sub)millimeter wavelength 
colors are rather red, and grains larger than ~1 mm are 
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needed to simultaneously fit the Herschel and 1.3 and 
3.3 mm fluxes. Instead, we find basically two ways to re- 
produce the observed far-IR Herschel colors: making the 
outer disk very small (by reducing the value of Rout) or 
very tenuous (by reducing the value of the surface den- 
sity profile, a). Both types of models contain little or 
no dust beyond ~40 AU, but each family has a different 
astrophysical implication: either T Cha has a small out- 
wardly truncated disk or material is accumulating at the 
inner edge of the gap. Resolved images are mandatory to 
break the degeneracy between a, and the inner and outer 
disk radii. The Atacama Large Millimeter/submillimeter 
Array (ALMA) will have both the resolution and sensi- 
tivity needed to directly measure the width of the disk as 
well as its surface density profile and distinguish between 
the two interpretations. The former result would encour- 
age even deeper searches for a wide-separation low-mass 
companion (most likely a Jovian planet). The latter sce- 
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TABLE 1 
T Cha Photometry data 



Wavelength 


Flux 


Error 


Telescope 


Reference^ 


(/xm) 


(mjy) 


(mJy) 






0.36 


6,80e+00 


10% 


ESO 50 cm 




0.44 


t: . w O 1 vJ ±. 


10% 


ESO 50 cm 




0.55 


1.12e+02 


10% 


ESO 50 cm 




0.64 


2 nip+n2 


10% 


ESO 50 cm 




0.79 




10% 


ESO 50 cm 




1.2 


A 2np+n2 

t: . w V3 1 V J ^ 


10% 


2MASS 


2 


1.6 


O . 1 W ^ 


10% 


2MASS 


2 


2.2 


l.lOe+03 


10% 


2MASS 


2 


3.6 


1.49e+03 


10% 


Spitzer 


3 


4.5 


1.32e+03 


10% 


Spitzer 


3 


5.8 


1.07e+e3 


10% 


Spitzer 


3 


8.0 


6.66e+02 


10% 


Spitzer 


3 


70 


5.06e+03 


15% 


Herschel 


4 


160 


2.97e+03 


15% 


Herschel 


4 


250 


1.78e+03 


20% 


Herschel 


4 


350 


1.06e+03 


20% 


Herschel 


4 


500 


6.60e+02 


20% 


Herschel 


4 


1300 


1.05e+02 


15% 


SEST 


5 


3300 


6.40e+00 


15% 


ATCA 


5 



Note. — «^ References: (1) Alcala et al. 1997; (2) Skrutskie et 
al. (2006); (3) Brown et al. (2007); (4) this work; (5) Lommen 
et al. (2007) 



TABLE 2 
Fixed Parameters from Oil 



Parameter Adopted value 



Stellar Tg// [K] 5400 

Stellar radius [R©] 1.3 

Stellar mass [Mq] 1.5 

Distance [pc] 100 

Inclination [deg] 60 

Grain size distribution slope, p —3.5 

Inner disk R^^ [AU] 0.13 

Inner disk Rout [AU] 0.17 

Inner disk M^^^^ [M©] 10-^ 

Inner disk amzn [/^m ] 0.1 

Inner disk amax [l^^] 10 

Outer disk amin [m"^ 1 0.01 



TABLE 3 

Model Grid and Fitted Outer Disk Parameters 



Parameter 


values sampled by grid 


best-fit model 


a = —1 modeP 


a = —3 modeP 


Surf, density prof, (a) 


-2.5,-1.75,-1.1,0.75,0.5 


-0.75 


-1 


-3 


Rin [AU] 


4, 5.5, 7.5, 11, 15, 18, 25, 30, 40, 50 


18 


10.5 


15 


Rout [AU] 


20, 40, 80, 150, 300 


20 


30 


300 


Mdust [Mo] 


(10, 3, 1, 0.3, 0.1, 0.01,0.001) X 10-4 


1.0x10-3 


2.0x10-^ 


1.8x10-^ 


Flaring index ((3) 


1.00, 1.07, 1.15 


1.15 


1.1 


1.1 


H50 (H at 50 AU) [AU] 


4, 6, 8 


4 


5 


6 


a^max [mm] 


0.01, 0.03, 0.1, 0.3, 1, 3, 10, 30 


3 


3 


3 



Note. — " The best-fit models fixing the value of a to —1 or —3 and fine-tuning the other parameters by hand 
guided by the results from the grid. 
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Fig. 1. — Top Panel: The full Spectral Energy Distribution of T Cha and a variety of models as labeled. The Oil model 
fits all previously available data but greatly overestimates the new 160, 250, 350, and 500 /im fluxes, illustrating the power of 
Herschel to completely rule out certain models. Lower Panel: The SED region sampled by Herschel. While the global best-fit 
model has a far-IR slope slightly steeper than observed, all disk models with a shallow surface density profile and a large outer 
radius, including the Oil model, result in far-IR colors that are much too red. The surface density profile of the outer disk and 
its inner and outer radii are degenerate parameters. Equally good fits can be obtained by either adopting a small outer radius 
(Rout = 30 AU, red dashed line) or a very steep surface density profile {a = —3, light blue dot-dashed line). 
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Fig. 2. — Top Rows: Bayesian probabilities for 6 of the 7 free parameters in our model for the entire grid (blue histograms) 
and the subset of models with Rout = 300 AU (red histograms). Due to space limitations, the /3 parameter is not included, but 
it presents a rather flat probability distribution. Models with large Rout clearly favor steep surface density profiles and larger 
H50 values. Bottom Row: 2-dimensional probability maps of the full grid. The highest probability regions are indicated in 
bright colors. Many models with R^n 18 AU and a steep surface density profile {a < —1.5) fit the data relatively well (left 
column). Similarly, a large number of narrow rings {Rin 18 AU and Rin < 50 AU) can reproduce the observed SED (middle 
column). Models with Rout ^ 40 AU are most favored by the data. Disks with shallow surface density profiles {a > —1.0) and 
large outer radius (Rout 150-300 AU) result in poor fits to the SED (right column). 



